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Abstract 

The main objective of this research is to 
analyze the effect of acoustic 
disturbances on a premixed flame and 
determine their role in the onset of 
combustion instabilities. Computations 
for the one-dimensional, unsteady 
combustion of a lean premixed methane- 
air mixture are performed. An acoustic 
excitation is introduced in the chamber 
and interacts with the flame front. Our 
results indicate that as the amplitude of 
the acoustic excitation is increased, the 
flame front position fluctuates rapidly. 
This phenomenon is even more intense 
when the frequency of the acoustic 
disturbance matches the fundamental 
frequency of the chamber. Our results 
suggest that the interactions between the 
flame and the acoustic excitation may 
result in flame extinguishment. In 
addition various passive control devices 
are tested and we found that the 
Helmholtz resonator with rounded inlet 
comers is the most efficient. 

Introduction 

Combustion instabilities arise 
in many practical devices, such as rocket 
engines, aircraft engines, gas turbines, 
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etc. Combustion instabilities are a 
concern in many engineering devices as 
they lead to malfunctioning of the device 
and a reduction in its efficiency. The 
sources of these instabilities are diverse 
and are dependent on fuel-oxygen 
mixture, chamber geometry, injection 
mechanism, etc. It is believed that, 
combustion instabilities occur due to 
sustained pressure fluctuations of 
acoustic nature in a chamber where 
combustion takes place. These 
fluctuations alter the flame, which in 
turn drives sound through unsteady heat 
release. 

It is well known that the 
resonance between a fluctuation in the 
heat release rate and an acoustic mode of 
the chamber is one of the main reasons 
for combustion instabilities. Prasad [1] 
investigated the effect of an external 
pressure perturbation on “premixed 
flames” using detailed chemistry. He 
solved the one-dimensional, time - 
dependent, nonlinear conservation 
equations. His results showed that in 
response to the externally imposed 
pressure waves, the flame moves to a 
location where heat release and heat loss 
are in balance. 

Peters and Ludford [2] 
concentrated on small variations of 
pressure and found that it was not 
affecting the flame significantly and 
hence concluded that substantially large 
variation in pressure (of the order of 
ambient pressure) might be needed to 
quench the flame. On the other hand, 
Ledder and Kapila [3] studied the 
response of a premixed laminar flame to 
pressure perturbations with emphasis on 
flame extinction. It was found that a very 
small variation in dynamic pressure has 
an influence on the burning rate, which 
could eventually lead to flame 
extinction. 
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Due to the importance of 
combustion instabilities for different 
combustors, many studies have 
addressed the modeling and control of 
combustion instabilities in lean premixed 
mixtures. Paxson [4] modeled 
combustion instabilities in an LPP 
combustor using a sectored one- 
dimensional code and studied the effect 
of an abrupt change of area in the 
combustor rig and how it expedites the 
phenomena of combustion instabilities 
using single-step kinetics. The author 
concludes that the position at which fuel 
is injected has a significant effect on 
acoustic instabilities. 

Frendi [5] concentrated on the 
coupling between the unsteady heat 
release rate from combustion and the 
acoustic energy of the combustion 
chamber. Furthermore, he also studied 
the effect of adding a Helmholtz 
resonator and/or baffles to the system. 
He found that the use of a Helmholtz 
resonator and/or baffles shifts the 
frequency of the chamber away from the 
unstable mode. 

In addition to the above studies, 
several other studies on different aspects 
of combustion instabilities have been 
carried out [6-14]. The basis for all the 
previous studies is due to Lord Rayleigh 
whose criteria states that the acoustic 
waves in the chamber will amplify if the 
acoustic pressure and heat release rate 
are “in phase”. 

In the present study, we 
concentrate on excitation of acoustic 
oscillations in combustion chambers, as 
it is an undesired phenomenon which 
results in inefficient burning and can 
lead to possible engine failure. However, 
accurate predictions of these types of 
acoustic instabilities are difficult to 
make as it primarily depends on the 
geometry of the combustor in question. 


The work performed describes the 
interaction of an acoustic excitation with 
a premixed laminar flame in a closed, 
one-dimensional chamber. Depending on 
the excitation frequency and amplitude, 
this interaction leads to the growth of the 
pressure fluctuations in the system. At 
substantially high amplitude, the growth 
becomes nonlinear and results in flame 
extinction. 

Mathematical Formulation 

The Mathematical model consist 
of the Navier-Stokes equations combined 
with a one-step irreversible exothermic 
chemical reaction of the form 


v f F +v q O ~>v p P (1) 


where “F” is the fuel, “O” the oxygen 
and “P” the products with v, the 

corresponding stoichiometric 

coefficients. All species are assumed to 
be in the gas phase. 

The full set of conservation 
equations in one dimension are: 
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In addition, the equation of state is given 
by 


P = P=T 
M 


( 6 ) 


The source terms in the fuel and oxygen 
species conservation equations are 
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The equations are solved in a 

nondimensional form using the 

following reference variables: 
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Method of Solution 


with p being the pressure, p the mixture 
density, T the temperature, R the 
universal gas constant and M the 
mixture molecular weight. 

The source terms in the energy equations 
are given by 
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The set of PDE’s are solved 
using an explicit finite difference 
scheme. The scheme, which is a 
generalization of the original 
MacCormack scheme obtained by 
Gottlieb and Turkel [15], is fourth order 
accurate in space and second order 
accurate in time. The numerical scheme 
applied to a one-dimensional equation of 
the form 
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consists of a predictor step given by 
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followed by a corrector step of the form 
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In the above equations, the 
subscript i denotes the spatial grid point 
and the superscript n denotes the time 
level. Superscript (* * ) corresponds to 
an intermediate time level (between n 
and {n + 1)). The fourth order accuracy 
is obtained by alternating the scheme 
given above with its symmetric variant. 
Furthermore, Jameson’s [16] artificial 
dissipation terms were explicitly added 
to the finite difference scheme, which 
captured sudden jumps in pressure and 
other variables near shock waves and 
stagnation points. Adding these terms 
helped prevent the appearance of the 
well-known Gibbs phenomenon in 
regions containing severe pressure 
gradients. Both the inflow and outflow 
boundaries were rigid, u = 0 . 


Results and Discussions 


Consider a lean mixture of 
equivalence ratio 0 = 0.662 in a closed 

chamber. Applying an external acoustic 
excitation at the inflow boundary of the 
form 


it = e sin( litft + <J>) (14) 


where, u is the non-dimensional 
velocity at the inflow boundary, / the 

frequency, e the amplitude and O the 
phase. 

The spark ignition parameters 
used for the present system are 
r li; =0.004w ( r lg =2.75 fls and 

2£* = 2.95* 10" J / m 2 s . The domain 

length is 0.2 m. Gas properties for a lean 
mixture are different. We obtain the 
properties for the lean side from Coffee 
and Kotlar [17] and values 
corresponding to <p = 0 .662 are given in 

Table 1. Two different cases are 
simulated based on changing frequency, 
a resonant case, /= 1940 Hz, and an 
off-resonance case, / = 1 730 Hz . 

1) Resonant Case 

By changing the amplitude, £ , three 
cases were studied, with the use of the 
resonant frequency / = 1940 Hz 

(i) £ = 0 - No Excitation 

(ii) £ = 0.03 - Low Excitation 

(iii) £ = 0.3 - High Excitation 

Several observations have been made 
from the following figures. Figures 1, 2, 
3 and show the temperature profiles at 
different times for £ = 0 , £ = 0.03 
and £ = 0.3, respectively. In Figure 1 we 
observe that the temperature distribution 
is symmetrical and the flame fronts are 
traveling towards both boundaries at 
equal speed and the temperature of the 
flame front is settling down to the 
adiabatic temperature of the lean 
mixture, which is T b = 2142. 2K . Unlike 
the profile in Figure 1, the profiles in 
Figures 2 and 3 are different. The 
acoustic excitation results in increasing 
the speed of the positive propagating 



flame and slowing the negative 
propagating flame. Pressure time Iraces 
for different excitation amplitudes are 
shown in Figure 4, which shows the 
growth of the pressure oscillations with 
time and by increasing the amplitude. 

2) Off-Resonance Case 

By changing the amplitude, £, three 
cases were simulated, with the use of 
off- resonance frequency of / = 1730/fe 

(i) £ = 0 - No Excitation 

(ii) £ = 1.2 - Low Excitation 

(iii) £ = 2.4 - High Excitation 

Figures 5, 6 and 7 show the 
profiles of temperature at different 
excitation amplitudes and it can be 
observed that the effect on the flame 
front location as the amplitude is 
increased is very significant. Flame 
fronts are observed to move with 
unequal speeds towards both boundaries. 
Increasing the excitation amplitude, the 
inflow flame front is moving with a low 
speed. 

Figures 8 compare the negative 
and positive propagating flame radius 
time histories for £ = 2.4 with £ = 0 . By 
increasing the amplitude (£), the flame 
location fluctuates significantly which 
could lead to possible extinction. 

Figure 9 shows the comparison 
of the time traces of the pressure at the 
center of the chamber at different times 
for the same amplitude (£=1.2) but 
different frequencies (/ = 17307/z and 
/= 19407fe). It is observed that 
pressure is growing much faster in the 
chamber when resonance frequency is 
imposed compare to the off-resonance 
frequency case. 


3) Passive Control 

Several passive control devices were 
tested numerically using a resonant 
acoustic disturbance alone, without 
combustion. The devices tested were, a 
baffle, a quaterwave and a Helmholtz 
resonator. In the later two cases, both 
smooth and sharp corners were 
considered. Figure 10 shows a 
comparison of the OASPL measured in 
the chamber the controlled and 
uncontrolled cases. It is clear from the 
figure that the Helmholtz resonator 
results in better control. Figure 1 1 shows 
the effect of rounding the various 
corners of the resonator. Our results 
show that rounding the inlet corners 
alone gives the best results. Figure 12 
shows that the Helmholtz resonator is 
less effective at high sound pressure 
levels. 


Concluding Remarks 

Based on the results presented in this 
paper, the following conclusions are 
made: 

1. The effect of the acoustic 
excitation is to displace the flame 
front. Depending upon the 
amplitude and frequency of the 
excitation, the flame front 
oscillations can be severe and 
may lead to flame extinction. 

2. When the frequency of the 
acoustic excitation is at 
resonance with the fundamental 
frequency of the chamber, we 
found that the combustion- 
acoustic interactions in the 
chamber resulted in rapid growth 
of the pressure oscillations. It can 
also be said that these nonlinear 
interaction between the heat 



release from the flame and the 
acoustic excitation may lead to a 
possible flame extinguishment. 

3. The Helmholtz resonator with 
rounded inlet corners is the best 
passive control device tested. 

Acknowledgements 

Support for this research was 
provided by a grant from NASA 
Marshall Space Flight Center, Grant 
number NCC8-200, with Tom Nesman 
as the technical monitor. 

References 


[1] Prasad, K., “Interaction of Pressure 
Perturbations with Premixed Flames,” 
Combustion and Flame, Vol. 97, 1994, 
pp. 173-200. 

[2] Peters, N. and Ludford, G.S.S., “The 
effect of pressure variations on premixed 
flames,” Combustion Science 
Technology, 34, 331. 

[3] Ledder, G., and Kapila, A.K., “The 
Response of Premixed Flames to 
Pressure Perturbations,” Combustion 
Science and Technology, 1991, Vol. 76, 
pp. 21-44. 

[4] Paxson, D., “A Sectored One- 
Dimensional Model for Simulating 
Combustion Instabilities in Premix 
Combustors,” Glenn Research Center, 
Cleveland, Ohio, 2000, AIAA-2000- 
0313. 

[5] Frendi, A., “On the role of Acoustic 
Coupling on Combustion Instabilities,” 
AIAA-2003-3181, Hilton Head, South 
Carolina, 2003. 


[6] Wu, X., Wang, M., and Moin, P., 
“Combustion instability due to the 
nonlinear interaction between sound and 
flame,” Center for Turbulence Research, 
Annual Research Briefs 2001. 

[7] Matalon, M., and Matkowsky, B.J., 
“Flames as gas dynamics 
discontinuities,” J. Fluid Mech. (1983), 
vol. 124, pp. 239-259. 

[8] Landau, L.D., “On the theory of slow 

combustion,” 1 944, Acta 

Physicochimica URSS 19, 77. 

[9] Umurhan, O.M., “Exploration of 
fundamental matters of acoustic 
instabilities in combustion chambers” 
Center for Turbulence Research, 1999. 

[10] Lee, D.H., and Lieuwen, T.C., 
“Premixed Flame Kinematics in a 
Longitudinal Acoustic Field”, AIAA-01- 
3851, Joint Propulsion Conference, 
2001. 

[11] Paschereit, C.O., Flohr, P., 
Gutmark, E., and Bockholts, M., 
“Numerically and Experimental 
Analysis of Acoustically Excited 
Combustion Instabilities in Gas 
Turbines” AIAA 2000-2076, 
Aeroacoustics Conference, June 2000. 

[12] Culick, F.E.C., “A Note on 
Rayleigh’s Criterion,” Combustion and 
Science Technology, Vol. 56, 1987, pp. 
159-166. 

[13] Nguyen, Q. V., “Measurements of 
Equivalence Ratio Fluctuations in a 
Lean Premixed Prevaporized (LPP) 
Combustor and its correlation to 
Combustion Instability” National 



Aeronautics and Space Administration, 
Cleveland, OH. 

[14] Harper, J., Johnson, C., Neumeier, 
Y., Lieuwen, T., and Zinn, B.T., 
“Experimental Investigation of the 
Nonlinear Flame Response to Flow 
Disturbance in a Gas Turbine 
Combustor,” AIAA-0 1-0486. 

[15] Gottlieb, D., and Turkel, E., 
“Dissipative Two-Four Methods for 
Time- Dependent Problems,” 1976, 
Mathematics of Computations, Vol. 30, 
No. 136, 703-723. 

[16] Jameson, A., Schmidt, W., and 
Turkel, E., “Numerical Solutions of the 
Euler Equations by Finite Volume 
Methods Using Runga-Kutta Time- 
Stepping Schemes,” 1981, AIAA- 81- 
1259. 

[17] Coffee, T.P., Kotlar, A.J., and 
Miller, M.S., “The Overall Reaction 
Concept in Premixed, Laminar, Steady 
State Flames I. Stoichiometries,” 
Combustion and Flame 54:155-169 
(1983). 



Figure 1: Temperature profiles in a “Closed 
Boundary System” as a function of time for 
a lean (0= 0.662) flame with no excitation 
(£ = 0), / = 1940 Hz 


Time (s) 



Figure 2: Temperature profiles in a “Closed 
Boundary System” as a function of time for 
a lean (0 = 0.662) flame with low excitation 
(£ = 0.03 ),/= 19407/z 



Figure 3: Temperature profiles in a 
“Closed Boundary System” as a function 


of time for a lean (0 = 0.662) flame with 
low excitation ( e = 0.3), / = 1940 Hz 



Figure 4: Time traces of pressure at the 
center of the chamber with different 
excitations fir a lean mixture in a closed 
chamber (a) e = 0 , (b) e = 0.03 , (c) e = 0.3 
and (d) £ = 1.2, /= \940Hz 



a lean (0=0.662) mixture with no 
excitation ( £ = 0 ), / = 1 730 Hz . 
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Figure 6: Temperature profiles in a “Closed 
Boundary System” as a function of time for 
a lean (0 = 0.662) mixture with low 
excitation ( e = 1 .2 ), / = 1730 Hz . 
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Figure 7: Temperature profiles in a “Closed 
Boundary System” as a function of time for 
a lean (0= 0.662) mixture with high 

excitation ( e = 2.4), / = 1 730 Hz . 




Tine (s) pressure level in the chamber obtained 

with the various passive devices. 

Figure 8: Flame trajectories as a function of 
time for (£ = 2.4) and 

(£ = O),0= 0.662,/ = 1730/fe. 




Figure 11: Effect of smoothing the 

resonator corners on the overall sound 
pressure level. 


Figure 9: Time traces of pressure at the 
center of the chamber with different 
frequencies for a lean mixture in closed 
chamber (a) / = 1 730 Hz , (b) / = 1 940 Hz 
for £ = 1.2, (j) = 0.662 . 
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Figure 12: Effect of a Helmholtz 
resonator on the overall sound pressure 
level of a combustion chamber at high 
noise levels. 


